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qn^ Abstract: We derive rates of convergence for limit theorems that reveal the intricate 

structure of the phase transitions in a mean-field version of the Blume-Emery-Griffith 

^ model. The theorems consist of scaling limits for the total spin. The model depends 

P^ on the inverse temperature /? and the interaction strength K. The rates of convergence 

r^ results are obtained as (/3, K) converges along appropriate sequences (/3„, Kn) to points 

"trt belonging to various subsets of the phase diagram which include a curve of second- 

^ order points and a tricritical point. We apply Stein's method for normal and non- 

i— ' normal approximation avoiding the use of transforms and supplying bounds, such 

,_^ as those of Berry-Esseen quality, on approximation error. We observe an additional 

K»- phase transition phenomenon in the sense that depending on how fast K^ and /3„ 

_^^ are converging to points in various subsets of the phase diagram, different rates of 

t~^ convergences to one and the same limiting distribution occur. 
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1. Introduction 

1.1 The Blume-Emery-Grifflths Model 

In 1971 Blume, Emery and Griffiths [2] introduced a mean field version of an important 
lattice spin model due to Blume and Capel. We refer to the mean field version as the BEG 
model. The BEG model is equivalent to the Blume-Capel model (see [T] and [3], |3j and |5j) on 
the complete graph on n vertices. One of the most outstanding features of the model is that it 
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is one of the few mean-field models that exhibits a continuous second-order phase transition, a 
discontinuous first-order phase transition and thus has a tricritical point, which separates the 
curves of the points belonging to the phase transitions. As a consequence this model is used to 
study many diverse systems, obviously including the one Blume, Emery and Griffiths devised 
it for. 

They showed that the model can be used to determine the phase diagram of He^-He^ mixtures 
using a simplification. In order to analyse this physical system the BEG model was introduced, 
which can also be used to explain the behavior of other physical systems such as microemulsions, 
semiconductor alloys or solid-liquid-gas systems, to name only a few. A variety of these appli- 
cations of the BEG model are discussed in jTHl Section 1]. Especially because the model keeps 
the intricate phase transition structure it continues to be of interest in statistical mechanics. 

Next we will give a mathematical definition of the BEG model and state some of the results 
known. 

Let /3 > and K > 0. As a configuration space for the model we will take all the sequences 
(cui, . . . ,Un) in { — 1,0, 1}". Ui denotes the spin on site i of the complete graph on n vertices. 
The Hamiltonian for the BEG model is defined by 

Hn,KH = Y.^^--[Y.^A (1.1) 

for each u G { — 1,0,1}". K > represents the interaction strength of the model. Given this 
Hamiltonian the probability of observing a subset A of {—1, 0, 1}" equals 

PlS,K,niA) = -77^ I exp(-l3H^,K)dPn. (1.2) 

Zl3,K,n denotes the normalisation constant and P„ is the product measure on { — 1, 0, 1}" having 
identical one- dimensional marginals p = |(5_i + (5o + (5+i). We will be interested in the behavior 
of the spin per site 

1 1 " 

-Sn{uj) := -y^ujj (1-3) 

under the distribution Pp,K,n as n — )■ oo. S'„ is called the total spin. The BEG model shares the 
feature - with for example the Curie- Weiss model - that the interaction terms in the Hamiltonian 
can be written as a quadratic function of the total spin. For this purpose we absorb the first 
non-interacting part of the Hamiltonian (1.1[) into the product measure. It is important to 



notice that in that case, in contrast to the Curie- Weiss model, the BEG model has a much 
more complicated product measure P^^p on { — 1,0, 1}" because of its /3-dependence. The one- 
dimensional marginals of Pn^p are given by 

pp{duj) = — - ■ exp{-f3ujf)p{dwj), 

where Z{(3) is equal to /exp(— /3w|)p(da;j) = ^^^ — . Thus one has (see (THl Section 3.1]) that 
the probability of observing a configuration u equals 



^/3,x,n(dcj) = -~ exp 



n^K r" 



SJu) 



\ n 



PnA^^) (1-4) 
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with normalization constant Zp^K,n = z(r)" • Hence one has reduced the BEG model to a 
Curie- Weiss-type model. 

We appeal to the theory of large deviations to define the set of (canonical) equilibrium 
macrostates. In order to state a large deviations principle (LDP) (for a definition see |Hl 
Section 1.2]) for the spin per site for the BEG model we need to define the cumulant generating 
function of pjs, which is given by 

f /l + e"^(e* -F e"*)\ 
C/3(t) = log J exp{tui)pp{dui) = log I ^ ^ ^^_p 1 . (1.5) 

Cramer's theorem ([SI Theorem 2.2.3]) states that, with respect to the product measure Pn,/3, 
the sequence (5'„/n)„ satisfies the LDP on [—1,1] with speed n and rate function 

Jp{z) := sup{te - C/3(t)}, (1.6) 

which is the Legrende-Fenchel transform of c^. Having the LDP for {Sn/n)n with respect to 
Pn,i3, the following theorem, taken from P^ Theorem 2.4], states the LDP for {Sn/n)n for 

Pl3,K,n- 

Theorem 1.1. For all /3 > and K > the following conclusion holds: With respect to Pi3^K,n, 
{Sn/n)n satisfies the LDP on [—1, 1] with speed n and rate function 

I^z) = Mz) - /3Kz' - mUMy) - ^Ky^}, 



with Jp{z) taken from (1.6). 



As a consequence only the points z G [—1, 1] satisfying Ip^K^z) = do not have an exponen- 
tially small probability of being observed. These points form the set of the so-called equilibrium 
macrostates, which is accordingly defined as 

Mp.K = [ze [-1,1] ■.Ip,k{z)=q}. (1.7) 

In [ISl Theorem 3.6, Theorem 3.8] it is proven that there exists a critical inverse temperature 
/3c = log 4 and, for /3 > 0, there exists a critical value Kc{fi) > characterising the phase- 
transition structure of the model: for /3 > and < K < Kc{(3), M/s^k consists of the unique 
pure phase 0, [ISl Theorem 3.6(a) and 3.8(a)]. For /3 > and K > Kc{(3), M.j3^k consists of 
two distinct, nonzero phases. For < /3 < /3c, as -ft' increases through -ft'c(/9), ■M.ii,K undergoes 
a continuous bifurcation, which corresponds to a second-order phase transition, [Ml Theorem 
3.6(b) (c)]. Here we have 

For /3 > /3c, as K increases through Kc{f3), Adp^x undergoes a discontinuous bifurcation, which 
corresponds to a first-order phase transition, [THl Theorem 3.8(c) (d)]. The point (/3c, Kc{f3c)) = 
(log 4, 3/ [2 log 4]) in the positive quadrant of the (3-K plane separates the second-oder phase 
transition from the first-order transition and is called the tricritical point. 

Based on the points that correspond to the transitions we define two different sets that will 
infiuence the form of our limiting density. The first set contains the single-phase region and it 
is defined by 
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The curve containing the second-order points is given by 

B = [{/3,K)eR^ : < (3 < /3c, K = K, 

K 




/3 



Figure 1.1: The sets A, B and C 

Furthermore we consider the singleton set 

C:=(/3„i^,(/3,)) = (log4,3/[21og4]), 

which separates the first- and second-order phase transition. See also [3, Section 1] for a nice 
summary. Figure 1.1 illustrates the sets. 

When studying the law of large numbers for the BEG model this complex phase-transition 
structure proves to be a determining factor. In [3 (2.1)] it was shown that for /3 > and 
< K < Kc{l3) the law of large numbers holds with Pii^K,n{Sn/n G dxj ^ 5q, as n — > oo. 
Hence, for sufficiently small interaction strength i^ > an analogue of the classical law of large 
numbers can be proven. If the interaction strength exceeds the critical value for K the law of 
large numbers breaks down. For /3 > and K > Kc{P) in [121 Theorem 3.6, 3.8] the existence 
of z{P, K) > was proven such that the following limit holds true: 

1 



Pl3,K,n[Sn/n G dxj => -{S^{f5,K) + S~ 



z{l3,K) 



see [3 (2.2)]. Because of the intricate phase transition structure there are also two limits for 
K = Kc{(5). Whereas for < /3 < /3c the law of large numbers holds, Pfi,Kc{p),n\Sn/n G dxj ^ 
5o, for (3 > f3c the limit is expressed by a measure supported at three points corresponding to 
the macrostates in A4j3^k'- 

where Aq and Ai are positive numbers satisfying Aq + 2Ai = 1 (for an explicit display see [H 
(4.4)]). These first hints of the intricacy of the phase-transition structure can also be seen for 
the limit theorems stated in Section |2l 

In Section [2] we will obtain limit theorems and rates of convergence for the rescaled total 
spin Sn/n^~'^ for appropriate choices of 7 G (0, 1/2]. In [7j 18 scaling limits and 18 moderate 
deviation principles for the total spin Sn were obtained as (/3, K) converges along appropriate 
sequences (/3„, Kn) to points belonging to the three separate classes: (1) the tricritical point 
C, (2) the curve B of second-order points, and (3) the single-phase region A lying under the 
curve. We obtain the 18 different scaling limits by an alternative proof (Stein's method) and 
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present rates of convergence in all 18 cases at the same time. Furthermore we observe that the 
complex structure of the phase transitions in the BEG model provides an additional insight, 
presenting that in 15 of 18 cases the rate of convergence differs within the same case. A fixed 
case - out of the 18 - is characterised by a fixed limiting distribution. We will observe that for 
any of the 15 different limiting densities the rate of convergence depends on the choice of the 
value of 7 and/or the choices of two more parameters Ai and A2 given as follows: We consider 
sequences (/3„, Kn) converging to (/3, K) taken from A, B 01 C along the sequences 

Pn = log (e/^^ - ^) , (1.9) 

K„ = K,{f3n) - 4", (1-10) 



where Ai > 0, A2 > 0, 6 ^ 0, A; ^ and K^{/3) defined in (|L8|) for /3 > 0. If {(3, K) is taken 
from a set A, B oi C the sequence (/3„, Kn) will determine the value of 7, since the sequence 
establishes which set influences the convergence towards (/3, K). Depending on the sign of h and 
k the sequences converge from regions having a different physical behavior. The mathematical 



explanation for the choices (1.9) and (1.10) will be clear later (the sequences are chosen so 
that certain terms in a Taylor expansion have appropriate behaviour), whereas the physical 
significance is not obvious. 

The three seeds from which the present paper grew are references [TB], [7j and [12]. In the 
first paper the phase transition structure of the BEG model is analysed. In the second paper 
limit theorems in the BEG model are proven and in the third paper, rates of convergence are 
obtained for limit theorems for the Curie- Weiss model when the inverse temperature converges 
to the critical inverse temperature in the model along an appropriate sequence (3n- These 
results generalise the limit theorems obtained in [T3] and [E]. The third paper developed 
Stein's method for exchangeable pairs for distributional approximations including the Gaussian 
distributions as well as non-Gaussian limit distributions and obtained convergence rates for the 
Curie- Weiss model (see also ^). Note that the fact that limit results are obtained as (3 converges 
along appropriate sequences /3„ is shared with a number of mean field models, including the 
Curie- Weiss models ([ID], [12]) and the Hopfield model of spin glasses and neutral networks 

([ig, [n]). 

1.2 The Function Gp^K and its Properties 

From now on we denote for a function / : M — ^ M the i-th derivative by Z^*-*. A crucial 
element for Stein's method is the function 

Gp^K^x) = PKx^ - Ci3{2(3Kx) (1.11) 



for X G M and its minima, cp denotes the cumulant generating function of pp given in (1.5[). 



The function Gjs^k plays a central role in nearly every aspect of the analysis of the BEG model, 
since it gives an alternative characterisation of the set of equilibrium macro states Aip^x- Apart 
from being helpful while developing Stein's method in the sequel, its usefulness is also certain 
for the study of, for example, the phase transitions, the LDP or moderate deviations for the 
total spin per site, see [7]. The fact that a wide variety of phenomena can be obtained via 
properties of a single function is an appealing feature which is shared with a number of other 
mean-field models including the Curie- Weiss model, the Curie- Weiss-Potts model [T7] and the 
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Hopfield model [12]. The next Lemma, proven in [TBI Proposition 3.4], draws a connection 



between the equihbrium macrostates defined in (1.7) and the minima of the function Gfi^K- 
Lemma 1.2. For each x eM. we define Gf^^xix) as in (1.11). Then for each /? > and K > 0, 

min< Jb(x) — BKx^ \ = tcvuiIGb k{x) k 

\x\<l^ PV / f- J ^gjj 1^ P,^v /J 



with Jpix) defined in (1.6). Additionally the global minimum points of Jib{x) — PKx"^ coincide 
with the global minimum points of G^^k o-nd thus 



■Mi3,K = |x G M : X minimizes Gi3^k{x) \. 



With the help of Lemma 1.2 the structure of these minima was discussed and proven in 
[m Theorem 3.6, 3.8]. For /3 > and < K < K^{f3), Mp,K = {0}, for < /3 < /3, 
and K = K^P), Mp^K = {0}. For < /3 < /3c and K > Kc0), there exists z{/3,K) such 
that Mi3^K = {±z{l3,K)}. For (3 > (3^ and K = Kd/S), there exists z{l3,Kc{l3)) such that 
MpK = {0,±z{f3,Kc{P))}. Finally for ^9 > /3c and /if > Kc{f3), there exists z('p,Kc{f3)) such 
th^tMp,K = {±z{P,K,{(3))}- 



A crucial element for the analysis of the model is the Taylor expansion of (1.11 ). For general 



(/3, K) we have, since Gjb^k is real analytic, that for the global minimum point 

G^^^'Uo) 
GM^) = Gp,k{0) + ^fvpa^^'- + O (x2'-+i) as X -> 0, 

(zrj. 

since G'L^(O) = 0. Here, r denotes the type of the global minimum point. In [TJ Theorem 4.2] 
the types were determined in the following theorem. 

Theorem 1.3. For all {f3,K) e AU B U C, Mpj^ = {0}. 

(1) For all (/3, K) G A, the global minimum point has type r = 1. 

(2) For all {(3, Kc{(3)) G B, the global minimum point has type r = 2. 

(3) For all C = {(3c, Kc{(3c)), the global minimum point has type r = 3. 



Remark 1.4. For the values of the parameters that were not considered in Theorem 1.3 the 
type of the global minimum points is 1, which is proven in [TBI Theorem 6.3]. 

This theorem will yield the Taylor expansion of Gfs^x if {P, K) is fixed and taken from one 
of the sets A, B or C. Next we deal with the sequences (1.9[) and (1.10) and the associated 



function G^„,k„. The following theorem yields three different forms of the Taylor expansion of 

GjB^^Kr,, see [3, Theorem 4.3]. 

Theorem 1.5. For 'j G M+ and for a positive bounded sequence (/3„, i^„) the following conclu- 
sions hold. Let for R> 0, \x\ < R. 

(1) For {(3n, Kn) — > {(3, K) G A the type of the minimum point G M.i3^k is r = 1 and there 
exists ^ = ^{x) G [— x,x] such that 

G,^,kM = '"f^ ' x' + A,^,kM^))x'. (1.12) 

The error terms A/^ k {^{x)) are uniformly bounded over n G N and x G {—R,R). We 
have Gfl,^ (0) = (^X^^)(^;:^^-^/^^^^) . 
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(2) For {f3n,Kn) — ?• {f3,Kc{f3)) G B the type of the minimum point G M.p^Kc{i3) is r = 2 
and there exists ^ = C,{x) G [— a;,a;] such that 

G,.,kAx) = ^£|^x^ + ^SkM.^ + B,„,,M^)).^ (1.13) 

The error terms Bp^^Kn{^{x)) are uniformly hounded over n G N and x G {—R,K). We 
have <„V„ (0) = ^^"^'Y."J;2)^'^'''^ ^^d cflj,^ (0) ^ for n -^ oo. 

(3) For (/3„, fi'n) — ;• (/3c, -ft'c(/9c)) G C" i/ie type of the minimum point G A^/3^,ii'^(/3^) z's r = 3 
and there exists ^ = ^{x) G [— x,x] such that 

^ / \ Gl3„,Kr,i^> 2 , ^/3„,X„(0) 4 , ^ftj.A'^lO) 6 , ^ /r-/ w 7 /i i,i\ 

G';3„,X„(X) = ^^ X' + -^^ X^ + -^^ X^ + Cp^^KSi{x)W- (1-14) 

r/ie error terms Cj^^^Kni^i^)) ^^^^ uniformly bounded over n G N anc? x G {—R,R). We 
have G^Ik^ (0) ^ and G^l^^ (0) ^ for n -^ oo. 



(4) Furthermore let us choose the sequence Kn as in (1.10), then in (2)-(3) we obtain 

GtKA^) = ^^C^' (1.15) 

with C^^ = ^'^„"s — )■ 2(3c- Moreover if we assume that (3n is chosen as in (1.9) we have 



Kc{l3n) 

in (2)- (3) that 



G'IkA^) = ^,C^' (1.16) 

with C(4) = fl^ ^ I see m (4.9), (4.10), (7.4) and (7.5)]. 

We next preview the contents of the present paper. In the next section, Section 2, we 
state the hmit theorems and the rates of convergence for the total spin per site. We are able 
to obtain 21 different limiting densities that result from the values of k and h defining the 
physically dissimilar regions of points. A, B and C, the sequence (/3„, K^) is converging from. 
Formulating the corresponding rates of convergence we will see that the 21 cases split into 42 



cases that result from the values of 7, Ai and A2 defined in (1.9) and (1.10) and thus depend 
on the speed at which {f3n,Kn) -^ {f3,K). The proofs of our Theorems will be presented in 
Section |4} They apply Stein's method, which shortly will be introduced in Section 3. 

2. Limit Theorems and rates of convergence 

We are prepared to state our results. Because of the intricate structure of the model we are 
able to find three different limit theorems if (/3, K) is assumed to be fixed and 18 scaling limits 



for the total spin as {l3,K) converges along the sequences (/3„, /^n) defined in (1.9) and (1.10) 
to points belonging to the sets A, B and C defined in Section 1. Let 

1^,:=-^ with 7 6 (0,1/2]. (2.1) 

2.1 Rates of convergence for fixed {/3,K) and {/3n, Kn) — ^ {/3,K) G A 

First of all we assume that (/3, K) is fixed. Limit theorems for the spin per site were first 
discussed in [7]. 
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Theorem 2.1. Let Wy be defined in (2.1) 



(1) For {P,K) & A we have 7 = 1/2. If Z^ is a random variable distributed according to 
the normal distribution A^(0,E1^^^,2) ^^ ^^ with expectation zero and variance "KiWh^), 



we have that 



sup 



P{Wi/2 <t)- p[Za <t)\<Li- n-^/^ (2.2) 



for some constant Li depending only on {(3,K). 
(2) For {l3,Kc{(3)) & B we have 7 = 1/4. If Zb is a random variable distributed according 
to the probability measure on M with density fsit) '■= L2-exp {—cf^) , with c = c{Wi/4) = 
{AE,{W^u))~^ , and L2 the appropriately chosen normalisation constant, we have that 



sup 



P{Wi/A < t) - p{Zb <t)\<L^- n~^/\ 



for some constant L3 depending only on {P,Kc{/3)). 
(3) For {f3c, Kc{f3c)) we have 7 = 1/6. If Zc is a random variable distributed according to 
the probability measure on M with density fcif) '■= La ■ exp {—dt^) with d = d{Wi/Q) = 
{6K[W^,q))~^ and L4 the appropriately chosen normalisation constant, we have that 



sup 

teR 



P{Wy, <t)- P[Zc <t)\<L,- n~"\ 
for some constant L5 depending only on {/3c, Kc{Pc)) ■ 



Remark 2.2. Theorem 2.1 shows that the rate of convergence is affected by the set containing 
(/3, K). For the BEG model, to the best of our knowledge, our results are the first ones, where 
the quality of approximation was estimated. In region A we found an optimal rate n~^/^ for the 
Kolmogorov distance known as a Berry-Esseen type result. We do not know whether the other 
rates are optimal. There is one case known in the literature, where the limiting density is of the 
form exp (— const. x^) and the rate of convergence is of order n"^/^: This is the rescaled total spin 
in the classical Curie- Weiss model at the critical temperature /3c, see [^ Theorem 3.8]. The 
technical reason is that the Taylor expansion of the corresponding function G^? is given by the 
Taylor expansion of tanh(-). In the same paper, generalisations of the Curie- Weiss model lead 
to bounds of order n~^^^'^''^ whenever the limiting density is of type exp(— const. x^'^). Note that 
Theorems 5.5, 6.1 (Case 1) and Theorem 7.1 (Case 1) in [7J follow from our result. We present 
a proof by Stein's method, avoiding the application of transforms. The limiting densities in [7j 

show moreover, that for {P,K) e A we have \imn^^E{W^/^) = (^^^^0))"^ - (2/3K)-^ for 
(/3, K)eBwe have lim,_oo 4E(W^i%) = (,,^2)'?4-e^) (see [3 (6.5)]). Finally for (/3e, K^i^c)) we 
have \im.n^oc 6IE(W^f/e) = ^ (see [3 Theorem 7.1]). The choices of our densities of the random 
variables Za, Zb and Zc will be explained in Section 3 in more detail. If the limiting density 
of a random variable is not known and especially the limiting moments are unknown, it is a 
remarkable advantage of applying Stein's method to be able to compare the distribution of a 
random variable W of interest with a distribution which inherits some moments of W, which 
characterise the limiting distribution. 

If (/3„, Kn) denotes a positive, bounded sequence converging to (/3, K) E A the situation of 
the scaling limits is as follows: 
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Theorem 2.3. Let (/3„, Kn) be an arbitrary positive, bounded sequence that converges to 



{l3,K) E A. Then we obtain 7 = 1/2 and the same result as in (2.2), Theorem 2.1. 



2.2 Six rates of convergence for {f3n,Kn) — ^ {f3,Kc{f3)) G B 

If {(3n,Kn) denotes a positive, bounded sequence converging to {(3,Kc{(3)) G B the situation 
of the scahng hmits gets more complicated. The form of the hmit theorem depends on the 
Taylor expansion of Gj^^^k^ i^ the neighbourhood of the global minimum point 0. This will 
become evident in the proof of the next theorem and is physically motivated in [7, Section 6] . 

Theorem 2.4. For fixed P G (0,/3c), let /3n be an arbitrary positive, bounded sequence that 



converges to /3 and Kn be the sequence defined in (1.10). Let W^ be defined in (2.1). Then by 
continuity of Kc{-) we have {f3n,Kn) — > {(3,Kc{f3)) G B. Assume that 

V = min(27 + A2 - 1, 47 - 1) = 0. 

and let 6{a,b) equal 1 if a = b and equal if a ^ b. Let Za^^a2 be a random variable distributed 
according to a densities of the form 

fa„a, (x) := C ■ exp [-{6{v, 27 + A2 - l)a,x^ + 5{v, 47 - l)a2X^)) , (2.3) 

for certain constants ai, 02 and C . 

(1) If '~f = 1/4 and A2 = 1/2 there exist explicit constants ai,a2 7^ (depending on (3n, Kn, 
KW^u and KWy^ presented explicitly in the proof) such that for a constant C 



sup 



P(Wy, < t 



PiZa„a,<t 



< Cn-^l\ 



(2) 7/27 = I-A2, 7G (1/4,1/2), A2 G (0,1/2), we take a^ = (2E(VK^2))-i > and 02 = 
and obtain with a constant C 



sup 



p(w^ < t 



P\ Zai,a2 ^ ^ 



<C 



n 



n 



1-47 



(3) // 7 = 1/4 and A2 > 1/2 we take ai 
constant C 



and a^ 



sup 



P(Wy^<t)-P(Za,,a2<t 



<c 




,7 G (1/4,1/3], 
,7 G [1/3,1/2). 

(4E(M/^JY4))~^ and obtain with a 

,A2G (1/2,3/4), 
,A2>3/4. 



Remark 2.5. We observe that the limit theorems depend on the value of 7 and on Kn through 
the speed A2. In the first case (1) the limit-density is exp(— aix^ — 02^^) and hence the case is 
influenced by regions A and B and the rate is n~i. The case corresponds to the critical speed 
A2 = 1/2. The coeflicient ai depends on the sign of A; 7^ and hence yields two different limit 
densities, whereas 02 > (both can be seen from the proof). 

In the second case, the limit-density is exp(— aio;^) with ai > 0. With A2 G (0, 1/2) it 
corresponds to a slow convergence of Kn to Kc{P). In this case only region A influences the 
form of the limiting density. But we consider converging in distribution to a normal distribution 
even though the non-classical scaling is given by n^~"' with 7 G (1/4, 1/2). Now we consider 
an additional phase transition phenomenon, since the rate of convergence depends on 7: the 
breakpoint is 7 = 1/3 and the more classical the scaling the better the rate of convergence. 
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In other words we see that for 7 G (1/4, 1/3] the influence of region A is getting weaker and 
weaker in the sense of a slower rate of convergence. 

In the third case the limit-density is exp(— 02^;'^) for any k E N. This case is linked to an 
influence of B. Here the speed A2 is most rapid. Anyhow, an additional phase transition 
phenomenon occurs: the rate of convergence depends on A2 and is getting best for A2 > 3/4. 
So if we choose a speed A2 > 3/4 we obtain the best rate n~^^^ or in other words: we can force 
the speed to obtain a rate of convergence which is optimal in comparison to the observation in 
Theorem [2l| For any A2 G (1/2,3/4) we have the rate n-(^2-i/2)_ 



Summarising we consider 6 different cases with 4 different limit densities (compare with 
Figure 4 on page 530 in [7j). In 5 cases (with the help of a certain speed up of Kn) we 
obtain the same rates of convergence as for fixed (/3, K), see Theorem 2.1 A phase transition 
phenomena persists in the case of a non-classical scaling. 

Remark 2.6. As shown in [71 Theorem 6.1] we focus only on the case of f = 0. If f > one is 
not able to obtain any limit theorem. For v < the authors in [7| obtain moderate deviation 
principles for the total spin per site. The fact that v is required to be zero becomes evident in 
Lemma 14.61 



2.3 Thirtytwo rates of convergence for (/3„,i^„) — > {/3c, Kc{/3c)) 

The last limit theorem of this section correspond to the case that the sequence (/3„, Kn) con- 
verges to the tricritical point {f3c, Kc{(3c)) = (log4, 3/[21og4]). Here we also need the sequence 
f3n — ;■ /3 taken from (1.9). 



Theorem 2.7. Let (3n and K^ he the sequences defined in ( 1.9[ ) and (1.10). Then {(3n, K^) — > 
{f3c, Kc{(3c)) ■ Let W^ be defined in (2.1). Given 7 G [1/6, 1/2] we assume that 



w = min(27 + A2 - 1, 47 + Ai - 1, 67 - 1) = 0. 

Then, if the random variable Zh-^^f,,^),^ is distributed according to the probability measure on 
with the density 

fbub2,b3ix) ■= C ■ exp (- ((5(w,27 + A2 - l)bix^ + 6{w,A-f + Ai - l)b2X^ + S{w,6-f - 1)633;^ 

(2.4) 

for certain constants &i, 62, &3 o-'^d C . 

(1) If'~i = 1/6 and Ai = 1/3, A2 = 2/3 there exists explicit constants 61, 62, &3 7^ (depend- 
ing on Pn, Kn, E,Wl,Q with i G {2, 4, 6} presented explicitly in the proof) such that for a 
constant C 



sup 



P{Wy, <t)- P(^6„fe,b3 <t)\< Cn-'/\ 



(2) //27 = 1 - A2, 7 G (1/4, 1/2), A2 G (0, 1/2) and Ai > we take hi = {2E{W^))-^ > 
and 62 = ^3 = and obtain for a constant C 

-^1-4^-Ai ^ ^ g (1/4^ l/'ilAi G (0, 1 - 37), 



sup 



P{W, <t)- P(Z,,,,,,,3 < t 



<C{n-^ 



n 



,7e(l/4,l/3],Ai>l-37, 
,7 G [1/3,1/2). 
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sup 



p{w^ < t) - p(z, 



biMM 



. <t 



<c 



(3) // 27 = 1 - A2, 7 G (1/6,1/4], A2 G [1/2,2/3) and Ai > 2A2 - 1 we take h = 
{2E,{W^))^^ > and 62 = ^3 = and obtain for a constant C 

'n^-^-^-^i , 7 G (1/6, 1/5], 1 - 47 < Ai < 27, 

n^-'^ ,7G(l/6,l/5],Ai>27, 

n^-^'^-^i , 7 G [1/5, 1/4], 1 - 47 < Ai < 1 - 37, 

n-T ,7G[l/5,l/4],Ai>l-37. 

(4) //7 = 1/6, Ai > 1/3 and A2 > 2/3 we take fe^ = 62 = and 63 = (6E(W^f/g))-i > 
and obtain for a constant C 

nV3-A, ^ ^^ ^ (1/3^ 1/2), A2 G (2/3, 5/6), 

Ai < A2 - 1/3, 
^2/3-A, ^ ^^ ^ (1/3^ 1/2), A2 G (2/3, 5/6), 

Ai > A2 - 1/3, 
^1/3-Ai ,Ai G (l/3,l/2),A2>5/6, 
^2/3-A2 ,Ai>l/2,A2G (2/3,5/6), 
n-i/6 ,Ai > 1/2,A2 > 5/6. 



sup 



P( W^l/6 < t) - P(^6„6„63 < t 



<C{ 



(5) 7/47 = 1 - Ai, 7 G (1/6, 1/4), Ai G (0, 1/3) and 2A2 > Ai + 1 u;e taA;e 61 = 63 
and 62 = (4]E(W^))^^ > and obtain for a constant C 

'n^-^-*-^^ ,7G (1/6,1/5),A2<47, 







sup 



p{w^ < t 



P Zi 



^bi,b2,i 



<t 



<c 



ni-67 ,7G (1/6,1/5),A2>47, 

ni-2T-^2 ,7G [l/5,l/4),A2<l-7 

-7 



n ' ,7G[l/5,l/4),A2>l-7. 

(6) If 'y = 1/6, Ai = 1/3 and A2 > 2/3 t(;e take hi = and 62,^3 7^ (depending on Pn, 
Kn, KWyQ with i G {4,6} presented explicitly in the proof) and obtain for a constant C 

^-{A2-2/3) ,A2G (2/3,5/6), 



sup 

teR 



P[Wy^<t)-P[Z,,^,,^,,<t 



<c 



n 



-1/6 



,A2>5/6. 

(7) If 'J = 1/6, Ai > 1/3 and A2 = 2/3 we take 62 = and 61,63 7^ (depending on /3„, 
Kn, KWl,Q with i G {2,6} presented explicitly in the proof) and obtain for a constant C 

^-(Ai-i/3) ,AiG (1/3,1/2), 



sup 



p{w, 



1/6 



< t 



P z 



■'biMM 



. <t 



<c 



n 



-1/6 



,Ai>l/2. 



7/47 = 1 - Ai, 7 G (1/6, 1/4), Ai G (0, 1/3) and 2A2 = Ai + 1 we take 63 = and 
61,62 7^ (depending on f3n, Kn, EWLq with i G {2,4} presented explicitly in the proof) 
and obtain for a constant C 



sup 



P{W^ <t)- P{Zi,„b,,bs < t 



<c 



n^-^T ,7 G (1/6, 1/5], 
n-^ ,7 G [1/5,1/4). 



Remark 2.8. We observe that the hmit theorems depend on the value of 7 and on /3„ and Kn 
through the speeds Ai, A2. In the first case (1) the hmit-density is exp(— 61x^ — 620;'' — 630;^) and 
hence the case is influenced by regions A, B and C and the rate is n^e . The case corresponds 



12 PETER EICHELSBACHER AND BASTIAN MARTSCHINK 

to the critical speeds Ai = 1/3 and A2 = 2/3. The coefficient 61 depends on the sign of A; 7^ 0, 
the coefficient 62 on the sign of 6 7^ and hence yields 4 different limit densities, whereas &3 > 
(both can be seen from the proof). 



The second case should be compared with the second case of Theorem 2.4: in addition to 



the conditions in Theorem 2.4 (2), we assume that /?„ converges with speed Ai to l3c- Again 
we consider convergence in distribution to a normal distribution even though the non-classical 
scaling is given by n^~'^ with 7 G (1/4, 1/2) and again the speed is n~'^ for any 7 G [1/3, 1/2), 
independent of Ai. But for 7 < 1/3 the rate of convergence is of order 77,1-^7- Ai g^^j hence 
slower than in region B. But if we speed up Pn in choosing Ai > 1 — 37 the rate of convergence 
is n~'' and hence faster than in region B. In total there are 3 subcases. 

The third case reads as follows. Under the same relation 27 = 1 — A2 as in case (2), with 
a speed up of Ai and A2 it is possible to observe normal convergence even for the scaling 
7 G (1/6, 1/4]. For 7 > 1/5 it is possible to obtain the rate n'"^ if we speed up (3n, for 7 < 1/5 
a speed up of /3„ implies the rate n^~^'^, which could have been expected in comparison to 



Theorem |2. 4 (2). Case (2) and (3) are linked to an influence of A. In total we have 4 subcases. 

The fourth and flfth case are linked to the limiting densities exp(— feax^) (influenced only 
by C) and exp(— 62a;'') (influenced only by region B), respectively. In case (4) a speed up 
of both Kn and (3n leads to the rate n~^^^, in case (5) the result is comparable with case 
(3): a certain speed up leads to rate n~'^ or n^~^'^ depending on 7. Interesting enough one 
obtains converging in distribution to exp(— 623^'^) even though the non-classical scaling is given 
by n^^'^ with 7 G (1/6, 1/4), which is comparable with case (2) and (3). Summarising we have 
9 subcases. 

The last three cases (6), (7) and (8) are linked to limiting densities exp(— 620:^ — b^x^) (influ- 
enced by B and C), exp(— feix^ — b^x^) (influenced by region A and C) and exp(— 6ix^ — b2x'^) 



(influenced by A and B). Cases (6) and (7) are comparable with case (3) in Theorem |2. 4 a 
certain speed up of Ai and A2, respectively, leads to the rate n~^^^. Finally case (8) is com- 



parable with case (2) in Theorem 2.4 with a non-classical scaling. In all three cases one of the 
nonzero parameters bi depends on the sign of fc 7^ and 6 7^ 0, respectively, and hence yields 
two different limit densities. Hence we have 4 cases each. 

Remark 2.9. As shown in [3 Theorem 7.1] we focus only on the case oi w = 0. If w > one 
is not able to obtain any limit theorem. For w < the authors in [7] obtain moderate deviation 
principles for the total spin per site. The fact that w is required to be zero becomes important 
in Lemma [4.61 

Summarising we consider 32 different cases with 13 different limit densities (compare with 
Table IV. on page 536 in [7]). In all but three cases (with the help of a certain speed up of /3„ 
and Kn, respectively) we obtain the same rates of convergence as for flxed (/3, K), see Theorem 



2.1[ A phase transition phenomena persists in three cases of non-classical scalings: cases (3), 



(5) and (8). 

Thus, combining the theorems above, we have 42 limit theorems depending on the values of 
(/3, K). Their proofs can be found in Section 4. 
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3. Stein's method 

Stein introduced in [H] the exchangeable pair approach. Given a random variable W, Stein's 
method is based on the construction of another variable W (some coupling) such that the pair 
(ly, Vr' ) is exchangeable, i.e. their joint distribution is symmetric. A theorem of Stein ( \21\ 
Theorem 1, Lecture III]) shows that a measure of proximity of W to normality may be provided 
in terms of the exchangeable pair, requiring W — VT to be sufficiently small. He assumed the 
condition 

E(I^'|iy) = (l-A)I^ 

for some < A < 1. Heuristically, this condition can be understood as a linear regression 
condition: if (PF, PF') were bivariate normal with correlation g, then E[iy|iy] = gW and the 
condition would be satisfied with X = 1 — g. Stein's approach has been successfully applied 
in many models, see e.g. [2T] or [22] and references therein. In [20], the range of application 
was extended by replacing the linear regression property by a weaker condition. We consider 
Stein's method by replacing the linear regression property by 

E{W'\W) = W + Xtlj{W)-R{W), (3.1) 

where iIj{x) depends on a continuous distribution under consideration and R{W) is a remainder 
term. Recently in [12] and [9] the exchangeable pair approach was extended to more absolutely 
continuous univariate distributions with a nice collection of new applications. 

Given two random variables X and Y defined on a common probability space, we denote the 
Kolmogorov distance of the distributions of X and Y by 

dK{X,Y) := snp \P{X < z) - P{Y < z)\. 

Let / = (a, b) be a real interval, where — oo<a<6<oo. A function is called regular if / is 
finite on I and, at any interior point of /, / possesses a right-hand limit and a left-hand limit. 
Further, / possesses a right-hand limit /(a+) at the point a and a left-hand limit f{b—) at the 
point b. Let us assume, that the regular density p satisfies the following condition: 

Assumption (D) Let p be a regular, strictly positive density on an interval / = [a,b]. 
Suppose p has a derivative p' that is regular on J, has only countably many sign changes, and 
is continuous at the sign changes. Suppose moreover that Jjp{x)\ log(p(x))| dx < oc and that 
ijj{x) := ^-7^ is regular. 

In [22] it is proved, that a random variable Z is distributed according to the density p if and 
only if K(f'{Z) + ipiZ) f{Z)) = f{b—)p{b—) — f{a+)p{a+) for a suitably chosen class J-" of 
functions /. The corresponding Stein identity is 

f\x)+ij{x)f{x) = h{x)-P{h), (3.2) 

where h is a, measurable function for which Jj\h{x)\p{x) dx < oo, P(x) := H^oViy) dy and 
P{h) := jj h{y) p{y) dy. The solution / := fh of this differential equation is given by 

,, , i:{Ky)-Ph)piy)dy 

/(-) = ^^ • (3-3) 



For the function h{x) := l{x<z}{x) let fz be the corresponding solution of (3.2). 
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Assumption (B) Let p be a density fulfilling Assumption (D) We assume that the solution 

f, of f',{x) + ^(x) f:,{x) = l{,j,<z}{x) - P{z) satisfies 

\Mx)\<d,, \f:{x)\<d, and \f'M-f:{y)\<d, 
and 



\m^)Ux)y\ 



,p'{x) 
' p{x) 



W))' 



< dA 



(3.4) 



for all real x and y, where di,d2, d^ and d^ are constants. 

We will apply the following results proved in [12]. Let pw be a probability density such 
that a random variable Z is distributed according to pw if and only if EfE[14^-?/'(14^)] f'{Z) + 

?/'(Z) f{Z)] = for a suitably chosen class of functions. 



Theorem 3.1. (see Theorem 2.5 in [12]^' Let p he a density fulfilling Assumption (D). Let 



iyV^W) he an exchangeahle pair of random variahles such that (3.1) holds with respect to p 
(ip = p' /p). If Zw is a random variahle distrihuted according to pw, we assume that the 



solutions f, of E[W^{W)] f'{x) + ^p{x) f{x) 
for any A > one has 



do 



^{x<z}\ 



P{z) fulfill Assumption (B). Then 



dKiW,Zw) < '^{Var{E[{W-W'y\W]y^\{d, + d2^Jnvn + ^A 



+ 



2X 

1( 
AV 4 



E(i?2 



d.A^ 



3A^ 



d?. 



+ —E{\i,{W)\) + ^E({W-W') 



^{\W-W'\>A} 



(3.5) 



Remark 3.2. In case the regression property (3.1) is fulfilled with i}) = p'/p, we expect a 



comparison of the distribution of W with Z distributed according to the regular Lebesgue- 
density p. Why do we introduce the modified density pvK? The reason was already discussed 
in [12]: If K^ is fulfilled, on obtains that E{W - W'f = -2XE[W^{W)] + 2E[WR{W)]. 



Comparing the distribution W with Z distributed according to p leads to a plug-in theorem 
(see [T^ Theorem 2.4]), where one has to estimate a term like 



E 



2A 



E[{W' -Wy\W] 



But with our observation E(l-^E[(l^'-iy)2|iy]) = l+E[Wij{W)]~lE[WR{W)]. Therefore 

the bounds in Theorem 2.4 in [12] are only useful, if -E[PF'?/;(iy)] is close to 1 and ^E[iyi?(iy)] 
is small. Alternatively, bounds can be obtained by comparing with a modified distribution that 
involves E[iy-?/'(iy)]. This leads to pw- Note that this is compatible with the quite general 
approach introduced in [H]. 

In the following corollary, we discuss the Kolmogorov-distance of the distribution of a random 
variable ly to a random variable distributed according to A(0,E(iy^)). 

Corollary 3.3. (see Corollary 2.10 in [12Jj Let a^ > and {W,W') he an exchangeahle pair 
of real-valued random variahles such that 



X 



E{W'\W) = fl - -^)W - R{W) 



(3.6) 
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for some random variable R{W) and with < A < 1. Assume that K{W'^) is finite and 
\W — W'\ < A for a constant A. Let Zw denote a random variable distributed according to 
N{0,E{W'^)). We obtain 



<y\^^ ,^._., ...,2„..i\Ai/2 J ^nW') {V2^ + A) \VE(i?2 



dKiW.Zw) < —iYi^Tmw'-WY\W])) +a2p^^ 4 - + 1.5A 



2AV V ^^ ' ' 'JJ ■ \ 4 X 



W -^ 1^^^^-^ + ^^^ j + ^' 1-5^ ^nm- (3.7) 

4. Proofs 

While developing an exchangeable pair {W^, W') and applying Stein's method for the BEG 
model we will be confronted with the conditional expectation of a single spin Wj and of the 
product of two single spins Ui Uj . Before proving the theorems we will collect some auxiliary 
results that will be needed in the sequel. The proofs will be quite elementary. Interesting 
enough the statements which follow will be the basis of our proofs. 



Lemma 4.1. Let uji G { — 1,0, 1}, Sn defined in (1.3) and S^^ := Sn — Wj- Then 

E [co,\Mk^,] = f^,KiSl^/n){l + Oil/n)) 

with 

„ , . 2e-^ smh{2f3Kx) .^^. 

^^'^^^^ '■= l + 2e-^cosH2PKxy ^^'^^ 

Proof. First of all we calculate the conditional probability of a given single spin 

P(3,K,ni{^i = n {(Wfc)fc^i}) 



Pl3,K,n{^i — t\{Uk) 



k^i) 



exp(-/3t2 + /^(t2 + 2tEwfe 



-1,0,1} ^ ^ kj^: 

exp{-l3e + ^{t' + 2tSl, 



ie{-i,o,i} ^ ^ ^^ 



(4.2) 



Thus we obtain 

IE[o"i|(afc)fc^i] 



e^2e-^sinh(2/3i^^ 



ie{-i,o,i} ^ ^ 



and with |t| < 1 it follows E[a.,\iak)k^i] < e^^^'/'^fp^KiSi/n) and E[ai\{ak)k^^] > 
e~'^/^^/"-fi3^j^(^Sl^/n). Hence the result is proved. D 



The next lemma will connect the function f^x defined in (4.1) and the function G^^k taken 



from (1.11). 
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Lemma 4.2. With the notions of Lemma \4.1\ we have 

fpAS'Jn) = Sl,/n-^Gf^(S\/n). 



Proof. A direct calculation and (|1.5|) yields 

d_ 
dx 

2(3K 



G^^'^j^ix) = ^{pKx'-c^i2PKx)) 



f 1 + 2e~P e ^^e"*""-" - e 

X — 



I3( 2l3Kx _ -2l3Kx 



(4.3) 



\^ I j^ (,-13^213 Kx j^^-2/BKx\ 1 + 26"^ 

= 23k(x 2e-/^sinh(2/3irx) \ 

^ \ l + 2e-('cosh{2f3Kx)J' 
which yields the result. D 

In order to get a bound on some variances we investigate the covariances for i ^ j. 
Lemma 4.3. For i ^ j , i,j G {1, . . . ,n}, we have 

EK^^|(c.,)z^{.,}] = /i(5;'Vn)(l + 0{l/n)) 

with 

, , , , , , 26-2/5(1 + cosh{ApKx)) 

fi{x) := fi,i3,K{x) :- 



1 + 4e-/5 cosh(2;gfs:x) + 2e-^l^{l + cosh{4:(3Kx)) 
and Sli^ := Yj ^t- Moreover we obtain 

nu^\{ui\^,]=f,{SlJn){l + 0{l/n)) 

with 

, , , , , , 2e-^ cosh{2f3Kx) 

/,(x) := /,,,,^(x) := i + 2e-/5,osh(2/3i^x) 

and S*^ = Sn — ^i- Also we have that < fi{x) < 1, i G {1, 2}, for all x E"^. 

Proof. First of all we take a look at the conditional probability of cjj and Uj given all the other 

spins. Pp^K,n[^i = t,UJj = s\{uJi)i^{ijy^ 

Pp,K,n ({Wi = t,^j = S}r\ {(Cj,),^{jj}}j 



Let 



exp(-/3(t2 + 8^) + ^ ((t + sf + 2(t + s)S\ 
E exp(-/3(m2 + r2) + ^('(m + r)2 + 2(m + r)S'A) ' 

m,re{— 1,0,1} 

A,.(a:) := exp(-/3(t2 + s^) + ^((t + sf + 2(t + s)x)), (4.4) 

Nt^s{x):=t\s^Dt,s{x)- (4-5) 
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Then we have that for the different values oi Ui = t and Uj = s 

t,sG{-l,0,l} 

We have 9 cases: For {ui,Uj) = (0,0) it is No^o{Si;^) = and £'o,o(5';'-') = 1. For 
{(0,1), (1,0)} we obtain A^o,i(>^^^) = and /^o,i(Sj^^') = e-'^+^^'^e^'^^""'^"'. If 
{(0, -1), (-1,0)} we obtain No^_i{Sf;^) = and Do,-i(5'^^') = e-'^+^^-'e-^/'^""'^"'. 
for {ui,uj^) = (1, 1) it is N^ i(5;'^) = e-2/3+4^"-^e4/^^""'^"', D^ ,{S'^^) = N^ i(5;'^'), for 
{(1,-1), (-1,1)} we have iVi,_i(^^^) = e'^^ = /^i,-i(^;'^). Finally for (u;„u;,) = (- 
holds iV_i,_i(5;'J) = e-^/^+^if'^-^e-^^-^""'^"' and D_i,_i(5;'^) = Ar_i,_i(^;'^). 

Using e^^-^""'^"' + e-^/^^""'^-' = 2cosh{4/3Kn-^Sii^) we have E[u;2u;2|(c^,),_,^.^^.j] 



(cJi,CJj) G 
Wi,Wj) G 

Moreover 

(wi,a;j) G 
-1,-1) it 



1 + 4e-/3+/3i^™-i cosh(2/3i^n-i5;'^) + 2e~^i^ + 2e-2/3+4/3^"-^ cosh(4/3i^n-i^^'^; 

2e-2/3 + 2e-2/3+^'^^""' cosh(4/3irn-i5;'J) 
1 + 4e-/3+/3i^n-i cosh(2/3i^n-i^;'^) + 2e-2/3 + 2e-2/3+4/3i^"-^ cosh(4;gi^n-i^;'^'; 

/i(5;'Vn)(l + 0(l^); 



Using the conditional probability of a single spin given all the other spins given in (4.2) we have 
exp(-/3 + ^(l - 25;)) + exp(-/3 + ^(l + 25, 



E 



UJi \{uJk)k^i 



ie{ 



E exp(-pP + ^(P + 2lSl 

-1,0,1} 



MSUn){l+0{l/n] 



D 



Lemma 4.4. Let 7 = 1/2 if (/3, K) e A, -f = 1/4 if (/3, K) e B and -f = 1/6 z/ (^, iT) = C. 
Then we have for i,j G {1, . . . ,n}, i 7^ j, 

Cov(cj2,a;|) = C(l/n'"^'^(^^'i)). 

where C denotes a constant. 

Proof. We have that 

Cov(w2,a;2) - ^ T^r, ,2, ,2,^ , ^ _ . il ^ T^r, ,2 



E 



EKMK^^)^^^^^^.^] -E E[a;2|^^^)^^^^^] E E[u^\{ui)i^{^y] 



Using Lemma 4.3 we obtain 



Gov (ulu^) = [Elf^iSi;^ /n)] - (E[/2(5;/n)])'j (l + 0(l/n) 

We observe that, for all i,j G {1, . . . ,n}, we have 

f2{SUn) = f2{Sli^/n) + U2{S\/n) - ^(S^'Vn)) = ^(S^'Vn) + 0{n~^). 

This follows since 1/2(3^) — /2(z/)| < 2e~'^| cosh(2/3i^x) — cosh.{2 (3 Ky)\ < c{l3,K)\x — y\, using 
Lipschitz-continuity of cosh(-) on a compact interval, where c{P,K) is a constant. Hence we 
obtain 

E[f,{S:i^/n)]-{E[f,{SUn)]y = E[f\{S'i;^/n)]-{E[f,{S:;^/n)^^^ 
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2_ Ae-^^ cosh\2/3Kx) 

~ 1 + 4e-/^ cosh{2/3Kx) + Ae-'^P cosh^{2PKx) ' 



We can see that 

(/2(X)) 

By applying the identity 2cosh^(x) = cosh(2x) + 1 we obtain 

if (x)? = 2e-^^{l + cosH4pKx)) ^ 

KJ2KX,, i + 4e-/3cosh(2/3irx) + 2e"2/3(i + cosh(4/3ii'a;)) ^^^'' 

Thus 

E[/i(S^V^)] - (E[/2(5;'Vn)])' = E[/|(5;'Vn)] - (e[/2(5;'V^)])' = ^MSli^ /n)]. 

Summarising we have 

CovK',c.|) = (N[h{S'^/n)] + nf2{Sli^ln)]0{l/n) + (9(1/^2)) (l + 0{l/n)). 
Since /2(0) = 0, Taylor expansion of /2 at leads to 

We note that 7 depends on the region of (;9, K): -f = 1/2 if {13, K) e A, -f = 1/4 if (/3, K) e B 



and 7 = 1/6 if {P,K) = C. Finally we obtain, by applying Lemma 4.6, that y[f2{S'!'^^ /n)] 
0{n-^^) and E[/2(^;'V^)] = 0{l) and therefore 

Cov(cj2,a;^) = C(l/n"""(^^'^)). 



D 



Remark 4.5. Note that a proof of Lemma 4.4 for parameters {/3n,Kn) converging to {/3,K) 



taken from one of the regions A, B and C follows exactly the lines of the proof of Lemma [4. 4| 
with j3 replaced by j3n and K replaced by Kn- For (/3„, Kn) the value of 7 depends on the region 
that the sequence is converging from. 

We can bound higher order moments as follows: 



Lemma 4.6. Let W^ be defined in (2.1). Then, for all positive bounded sequences (PnjKn), 
7 G (0, 1/2] and all I &N we obtain 

E[W{] < const. (l). 

Proof. The proof is based on the Hubbard-Stratonovich transformation that is used for example 
in [3 Lemma 4.1] to derive the central limit theorem for the total spin per site. The situation 
of fixed {(3,K) is included in the study of sequences {l3n,Kn) converging to {(3,K). Let Yn 
be a sequence of A^ (0, (2/3„i^„)~^) random variables independent of all other random variables 
involved. According to [3 Lemma 4.1] we have 

/ Y \ 1 

where Z denotes the normalisation. Obviously this transformation does not change the finite- 
ness of any of the moments of W.y. In order to use a Taylor expansion we have to differentiate 
between the regions A, B and C. If the sequence converges to a point in A we have 7 = 1/2 



and by (1.12) the density with respect to the Lebesgue measure is given by const. exp(—?/^) (up 
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to negligible terms). Next we turn to the set B. We consider an arbitrary positive, bounded 



sequence converging to (3 and Kn given by (|1.10|). With (|1.13|) we obtain 

1 fcC?) 



nG 



I3„,K, 



{y/rP 



n 



27+A2-1 



^" -y' + 



1 



^47-1 24 n^T"! /^n.-'^r, 



(O2/' 



where C^^-* — )■ 2/3^, see (1.15). In order to obtain a density that is given with respect to the 

-y^ — y*) (up to negligible 



Lebesgue measure by const.exp(— y^), const. exp(—?/^) or const. exp( 



terms) we note that f = is required with v defined in Theorem 2.4, It remains to discuss the 
set C. According to (1.14) we have 



nGf- 



,K, 



{y/n') 



^27+A2-l 2 



y^ + 



1 



1 G 



I3„,K„ 



^67-1 



6! 



(0) 6^ 

— y + 



^47+Ai-i 24 
1 



C^-^n ^4 



77-1 



Cf 



n 



n 5-"-i 



XOi, 



— )■ |, see (|1.16|). In order to obtain a density that is given with respect to the Lebesgue 



with C^4) 

measure by const. exp(— (?(?/)) (up to negligible terms), where G{y) is a linear combination of 
the terms y"^, y^ and y^, we note that w = is required with w defined in Theorem 2/7 In 
each of these cases discussed above a measure with the density stated there has moments of 
any finite order. D 



We now consider the construction of an exchangeable pair {W^^ W^') in our model for W^ = 
-f^ = T^ Y17=i ^i' proving an approximate regression property. We produce a spin collection 
u' = (Wi)i>i via a Gibbs sampling piocedvue: select a coordinate, say i, at random and replace Ui 
by u)'^ drawn from the conditional distribution of the z'th coordinate given (cjj)jyj, independently 
from Wj. Let / be a random variable taking values l,2,...,n with equal probability, and 
independent of all other random variables. Consider 



7 



W^ - — ^ + -^ 

T ^1-7 „l-7 



n^ 



n^ 



1 ^ Uj'j 



n^-T 



Hence (1^7, W^) is an exchangeable pair and W^ — W!^ 
obtain 



iJj—LJj 

ni-7 ■ 



For F := cr(a;i, . . . , a;„) we 



E[iy^ - iy;|j^] 

With Lemma 14.11 and Lemma 14.21 we have 



11" 1 

2 = 1 



1 1 " 

2=1 



¥.[W^-W'\T] = -W^ 



n 



11" 



^^G^^^j,iSUn)){l + Oil/n: 



Using 



1 1 v^ ^ 



1 Sr, 



nn 



1-7 



1 ^n 

n^ 77,1-7 



-vf; 



n 



n 



,H^7, 



we obtain 



E[W^ 



w^;i-^] 



11" 1 

^,E2^4!k^»(i+c?(i/-))+0K/-^). 



(4.6) 



(4.7) 
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Alternatively with f[s,K{Sn/n) = fi3,K{Sn/n) + fp,K{Sl^/n) - ffs,K{Sn/n) we see 

1 1 



E[W^ - W'\T] 



ni-7 2I3K 



G^^!K{Sn/n){l + 0{l/n)) + 0(w,/n') + R^^k,, (4. 



with 



1 1 



Rp,Kn := -Tz:;-T.{hK{SUn) - f^,K{Sn/n)){l + Oil/n] 



1=1 



(4.9) 



Before proving the theorems we fix an easy but very useful bound on Ri3^k,j'- 
Lemma 4.7. There is a constant C depending only on /3 and K such that 



Proof. The denominator of f^^K given in (4.1) is larger than 1. Hence for any x,y G [—1, 1] we 
obtain 

\fpA^) - fpAy)\ < 2e-^| smh{2{3Kx) - smh{2/3Ky)\ 

+4e-2^| sinh(2;girx) cosh{2/3Ky) - smh{2PKy) cosh(2;gis:x)| 
< c{(3, K)\x -y\+ 4e-2^| sinh(2/3ii'(x - y))l 

using Lipschitz-continuity of sinh(-) on a compact interval and the hyperbolic Pythagoras, 
where c(/3, K) is a constant. It follows that 

IfMSnM - fMSJn)\ < -ciP,K)+Ae-'^{2PK/n) + Oin~'). 



D 



Proof of Theorem, 2.1. We will only prove part (1) and (3) of the theorem. The proof of (2) 
follows the lines of the proof of part (3). In each of the cases the exchangeable pair is constructed 



via a Gibbs sampling, see (4.7). We start with the proof of part (1). In order to be able to 
apply Corollary 3.3 we need the linear regression condition given by ( |3.6 ). Since {(3,K) E A 
we have 7 = 1/2. With the Taylor expansion of G'/?,^' given in (1.12) we have 



n 



G^^^ASJn) = ^Ggko) + OiiSlJnf), 



and therefore applying (4.6) we obtain with ( |4.7 ) 

1 1 



/2 - yyi/2\'^\ = n2BK^'^'^^^^^^^^ ^ ^^'^ 



nWi/2-W[,^\7] 



A 



a' 



Wi/2 + -Rl/2 



with A = i, a^ = ^1^ and 



Ry2 = Oi^Wl^ + 0(-^Wy2 



Hence we have (3.6) and can apply Corollary 3.3 In Remark 4.8 we will comment that a in the 



linear regres sion we found is not automatically the variance of the limiting normal distribution. 
By Lemma 4.6 we know that E(iyf/2) is bounded and therefore X'^JeJ^R^^^) = C(n~^/^). 
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Therefore the last three terms in (3.7) can be bounded by a constant (depending on /3 and K) 



times n ^/^. Next we have to consider the variance of 



E 



1 " 

(iy{/2 - 1^1/2) W1/2] = -^Y.^[{uj[-UJ,)^\Wi/2 



n^ 



n' 



It 



i=l 



- 2ujM 



(^'i\Wi/2 



=: A1 + A2 + A3. (4.10) 

To get an estimate of the variance of this expression we will bound the variances of the Ai and 
start by taking a look at the variance of Ai. 

Applying Lemma 4.4 with 7 = 1/2 we have Cov{ujf,ujf) = 0{n^^). This leads to the bound 



V [^i] = O {n~^). A conditional version of Jensen's inequality yields V [^2] < V [Ai]. Thus the 
variance of A2 has the same order as the variance of Ai. Furthermore Lemma [4. 1[ Lemma 4.2 
and Lemma 4^ yield 



\A. 



n^ 



i=l 



/ 1/2 \ / -^ \ 

V ^2 y ^ri'^y 



As a result of Lemma 'LQ the variance of ^3 can be bounded by a constant times n~^. Sum- 
marising these estimations the variance of E (W[,2 — Wi/2y\Wi/2 can be bounded by 9 times 

the maximum of the variances of the terms Ai, A2 and ^3, which is a constant times n~^. 
Thus, finally 

,2 



a 
2A 



V 



E 



{Wi/,-Wy2y\Wi/2\] =0{n-'/') 



which completes the proof for the region A. 

Next we turn to the region C, part (3) of Theorem 2.1, As has been said before, the proof 
for region B follows the steps of the next lines except that slight changes regarding the Taylor 



expansion have to be made. In order to apply Theorem 3.1 we have to show that the linear 



regression condition (3.1) is fulfilled. Applying the Taylor expansion of G/3^^Xc{/3c) ^^ (1-14), with 
d^Sl) and (|49| we obtain 



E[Wi/e-WU\r 



1 '^B.,K, 



Ci-^^c\Hc 



n 



^/^5\2f3,KM 






Thus 



nwy,-w[jr 



-AV^(H^l/6)+/?l/6 



with \ = n ^/•^ and ip{x) 
have used Lemma [4. 6[ 



G' 



(6) 

5!2/3c-ft:c(/3c) V6 



,(0) 



R 



/3c,Xc(/3c),l/6 



0(n ^^/^j, where we 



Note that 



i(w. 



EUyi/6t/'(Wi/6) 



m^) 



. Applying Theorem 



3.1 



we will compare the distribution 



This density as well as the density p with ip = p' /p fulfill assumption (B), see [T^ Lemma 2.2] 



of Wi/Q with a distribution with Lebesgue-probability density proportional to expf 
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Given the linear regression condition we are able to compute the bound given in Th eorem |3.1 

1 



Since \Wi/q — W[,q\ < ^^ we set A := ^^ and obtain with the help of Lemma 



4.6 



'd.A^ 



3A 



E 



HW: 



l/6j 



0{n 



-5/6 



using Lemma 4.6, we obtain that 



The last term in (3.5) is zero. From Lemma 



4.7 



we know that \R 



di + d2JE{W' 



n- 



E[i?2 



/6J 



i3c,K^{l3^), l/6\ 



O {n-"^) . 



0{n"^^l^) and 



Additionally 



E 



(W^i/6 - W^i'/e) Wi/6 



1 



2/3 



n 



{Ai + A2 + A^ 



with the AiS defined in (4.10). In order to be able to estimate the variance of this expression 

1/6 (since {13, K) e C) 
A conditional version 



again we will bound the variances of the Aj's. With Lemma 4.4 for 7 
we have Cov{uf,uj'j) = C(n~^/^) and therefore Y[Ai] 



of Jensen's inequality yields V [A2] < V[yli]. Thus V 
Furthermore Lemma |4.1[ Lemma |4.2| and Lemma 46 yie 

^8/3 



n-2/3Ai 

d 



V 



n 



-2 



2^2/3 



1 " 



U,K{SJn)(l + 0{l/n 



O 



8/3 



n 



+ 



l^M 



1 



n^ 



0[n 



As a result of Lemma 4.6 the variance of ^3 can be bounded by a constant times n ^. Sum- 

3 bounded by 9 times the max 
which is a constant times n~^. Thus, 



marising the variance of E (Wi/e — W^i/e) Wi/6 can be bounded by 9 times the maximum of 



the variances of the terms ^73^41, 



finally 



^^2 and 



:^273^3) 



2A 



V 



E 



{Wi/,-Wi/,)^\Wi/e]] =0{n 



-1/3 



which completes the proof for the region C. 



D 



Remark 4.8. Note that a"^ in the proof of part (1) is not the limiting variance in [Tj Theorem 
5.5]. The variance is {Cf^Ki^))-^ - {2/3K)~^ with Gf^j^{0) = 2/3K{l -2/3Kcf\o)). Interesting 
enough in the classical Curie- Weiss model, the prefactor a^ in the regression identity coincides 
with the limiting variance, see the proof of Theorem 3.7 in [12]: Here the limiting variance is 
(1 - /3)-i and the prefactor is (3/Gf{0) = (3/{(3{l - /3cf\o)). In the Curie-Weiss model one 
has cf{0) = 1 and hence (3/Gf\o) = {1 - (3y\ 



Proof of Theorem 2.3. Since zero is a unique minimum for the whole set A the proof requires 
exactly the same steps as the proof of part a) of Theorem 2.1 D 



Now we turn to the theorems involving the sequence (/?„, K^) that converges to (/3, K). 



Proof of Theorem 2^. Our goal is to apply either Corollary |3.3| or Theorem 3J_, depending on 
whether there is a Gaussian or a non-Gaussian limit. Given W^ again we construct a coupling 
W via Gibbs sampling such that {W^, W') is exchangeable. This will be used in order to 
calculate A and R to get the linear regression condition, which is, also due to the form of the 
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plug in the Taylor expansion ( |1.13[ ) to obtain 
with 



limit density, either taken from (3.1 ) or (3.6). Let J-" := cr(cji, . . . ,cj„). We start with (4.8) and 

1 



w',\:f] 



.(2) 



1. 



^Pn-^r 



gZk. (0) -w, + 



.(4) 



n 



mPnKr. 



gZkA^) zj^^w' + R' 



n 



«-'-'5(§) + '5(;S) + «^.-~ 



7' 



where Ri3„,K„,-y is defined in (4.9). We can use (1.15) to obtain 

A: 1 „. 1 



E[w^ - wm 



-w^ + 



GfU (0) 



-w: 



R^. 



(4.11) 



(4.12) 



Kc{Pn) nl+^2 -^ mPnKn ^"'''" nl+27 ^ 

Proof of part (1): Depending on the infiuence of regions A and B there are different expres- 



sions for A due to an application of either Theorem 3.1 or Corollary 3.3 We note that the 



choice A2 = 27 seem to be necessary to get the expressions prior to W^ and PV^ of the same 
order. We do not expect the choice 7 = 1/2 or 7 = 1/6: remember that the first summand on 



the right hand side of (4.14) lead to a Gaussian limit in the case of choosing the scaling 1/n 



whereas the second summand lead to a limiting density in the case of the scaling l/rt'/^. Hence 
l/m? should be expected to be overdesigned whereas 5/3 = 1 + 27 gives 7 = 1/3, which is at 
least a non-classical scaling. Hence we consider 7 = 1/4 and A2 = 1/2 and expect an infiuence 
of both regions A and B (and it is known from (TJ Theorem 6.1] that this is the right choice). 
So we end up with 



with A = ^ and V'(x) = -]^;p;;)a; - 
where we used Lemma [4.61 



-A^(H^l/4)+i?l/4 



1 



mPnKn 



GZk^ (0) ^' and R 



1/4 



R 



i3„,Kn,l/4: 



0{n-y^), 



Note that —. — "^^ ^ = '^^^^^^^^ with explicit formulas for ci = ci(/3„, Kn), C2 = C2(/3„, Kn) 



and C3 = C3(/3„, Km'E.{W^u),'E{W^u))- Applying Theorem 3.1 we will compare the distribution 






of Wi/4 with a distribution with Lebesgue-probability density proportional to exp 

^^^r^j- This density as well as the density p with ijj = p' /p fulfil assumption (B), see 
Lemma 2 .2] G iven the linear regression condition we are able to compute the bound given in 
Since \Wi/4 — W[u\ < ;^ we set A := ^^ and obtain with the help of Lemma 

i\ Q /I , , 



Theorem 



3.1 



1 /^^ 
A I 4 



3^^ 

+ — E 
2 



HW, 



l/4y 



The last term in (3.5) is zero. From Lemma 4.7 we know that \R/3„^k„,i/4:\ = 0{n '^^^) and using 
Lemma 14.61 we obtain that 



rfi + d2^j¥.{W^) + -A ] A-\/E[i?2 



/4J 



O {n-^/') 



Additionally 



E 



{Wy^-W[,:)''\Wy, 



n 



1/2 



{Ai + A2 + A3 
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with the Aj's defined in (4.10 ) except that the expectation is now taken for the measure Pn,/3„,x„- 
In order to be able to estimate the variance of this expression again we will bound the variances 
Keeping Remark 4.5 in mind we can apply Lemma 4.4 for 7 = 1/4 and get 

: 0(n-^) 



of the y4j's 

)2 , ,2 



Gov (cur, w; 



0{n ) and therefore V [Ai 



A conditional version of Jensen's 



inequality yields Y [A2] < Y[Ai]. Thus V n-^/^Ai 



V 



n 



Lemma 4.1| and Lemma |4.2| yield 



1 1 

2ni72 



I A, 



1 






-1/2^, 



0(^ 



0{n ^). Furthermore 



+ 



n^ 



As a result of Lemma -LG the variance of A^ can be bounded by a constant times n ^. Sum- 
marising the variance of E {W1/4, — W[u)'^\Wi/4 can be bounded by 9 times the maximum of 

the variances of the terms ^72^1, ^72^2 and ^72^3, which is a constant times n""'. Thus, 
finally 

^2 



2A 



V 



E 



{Wy^-Wi,,)^\Wy 



O 



n 



-1/2 



which completes the proof of part (1) where both regions A and B infiuence the limit distribu- 
tion. 

Proof of part (2): If the infiuence is from region A the term involving W^ has to be of smaller 
order than the term of W^. Hence, we note that the condition A2 < 27 has to be fulfilled. In [3 
Theorem 6.1] it is proved that the only interesting choice for 7 and A2 is to take A2 G (0, 1/2), 
7 G (1/4, 1/2) and 1 — 27 = A2. Let us discuss this case in detail. We consider the linear 
regression condition 

k 1 



with 



E[W^ 



R^ 



K\^] 



S\2PnKn 



KciP. 



n 



l+A: 



-Wv + R^ 



G 



(4) 



^^)zx^,^' + R^ 



n^ 



and R^ given in (4.11). Hence A 



11+^2 



and ip{x) 



KciM 



X 



-x/cr^ and we compare 



the distribution of W^ with a N{Q,¥.iyV^)) distribution. Since \W^ - W'^\ < 
A : = 

A^/\ 



1 

= n 



and obtain with the help of Lemma 4.6 that a'^l.hAJ¥,{W^) = 0{ri 



-2+37 



n' 



— we set 
M. With 



^ the second last term in (3.7) has the same order. From Lemma 4.7 



we know that \Rp„,Kr,n\ 
Summarising we have 



C»(n^-2) and therefore 0{R^ 



n 



-mm(2-7,i+27)^ uslug Lemma 



4.6 



A" 



nm 



On 



- min(7,47— 1) 



As we can see, the third order term of the Taylor expansion of Gf^^^K„ now infiuences the order 
of the remainder. We have 

1 



E 



(w^-w'Yiw. 



1-27 



n 



{A, + A2 + A3 



with the Ai's defined in (4.10 ). We can apply Lemma 4.4 to obtain Cov(a;f , cj| 



and therefore as seen before V 
Furthermore Lemma 4.1 



and Lemma 4.2 



V 



n 



-(1-27)4 



(-)(^^-mm(47,l)j 
(-)(^^-min(4,5-47)) ^ ^(ri-S+^T), 



yield 1^1^31 = 0{^) + 0{^). As a result 
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variance of E 



of Lemma |4.6| the variance of A^ can be bounded by a constant times n ^. Summarising the 

can be bounded by a constant times n"^^'^'^. Thus, finally 



{w^-w:^y\w-y 



a 
2A 



V 



E 



(W, 



w;y\w^ 



0(n 



2-27-5/2+27 



O (n-V2) . 



The case A2 G (0, 1/2) corresponds to the slowest convergence of i^^ — )■ Kc{f3) with (/3, Kc{f3)) G 
B, in which only A influences the form of the limiting distribution, which has a Gaussian 
density even though a non-classical scaling is given by n^"'''. We obtain an additional and 
remarkable phenomenon: for any 7 G (1/4, 1/3] the rate of convergence is l/n^'^~^ whereas for 
all 7 G [1/3, 1/2) we obtain the rate 1/n'^. 

Proof of part (3): Finally we consider the case which corresponds to the largest value of A2, 
namely A2 > 27. We take A2 > 1/2 and 7 = 1/4 and thus the most rapid convergence of 
Kn — )■ Kc{(3). Now we end up with 



E[i^i/4 - w[jr 



with A 



1 

„3/2 



and ip{^) 



1 (^(4) 



0{n~'^^^), where we used Lemma 
obtain with the help of Lemma 
and third term of order O (n~^^^ 



4.6 



4.6 



/4K J - -X^{Wi/4) + Ri/i 

(0) x^ and R1/4 = R^„,k„,i/4 + 7^ 
Again \Wi/4 - ^^1/4! < 
that the first summand in (3.5) is of order O (n~^^'^ 

we know that \Ri3„,k„,i/4: 
Using Lemma 



^, we set A := ^ and 



4.7 



From Lemma 
The second summand of R1/4, is of order 0{n^^^^^] 

mm(i/4,A2-i/2)^^ This is also an interesting phase 



( 



4.6 



0(n-7/4). 
we obtain that 



IA]X'^ 



E[Rl 



Oh 



transition: for any A2 G (1/2,3/4) we obtain a slow rate of convergence n~^'^^^^^'^\ but when 
Kn converges more rapid in the sense of A2 > 3/4, we obtain the rate n~^^^. The proof via 
Stein's method gives the information that in case (3), we have to assume A2 > 1/2. D 



Proof of Theorem 2/I_. Again our goal is to apply either Corollary |3.3| or Theorem 3J^, depend- 
ing on whether there is a Gaussian or a non-Gaussian limit. Given W^ again we construct a 



coupling W via Gihhs sampling such that (W^, PF') is exchangeable. We start with (4.8) and 
plug in the Taylor expansion (1.14) to obtain 

1 



E[W^ 



w',\:f] 



G 



(2) 



+ 



WnKn 
1 



(0) 



"'^" V-^ -W, + ^y^^-^ 



G 



(4) 



,Kr 



(0) 



n 



1+27 "^'7 



.(6) 



GZk,. (0) 



1+47 



W^ + R^ 



n 



with 



R^ 









1+57 



n 



+ R 



■l3„,K„,-f, 



where RB„,Knn ^^ defined in (4.9). We can use (1.15) and (1.16) to obtain 



¥.[w^ - wm 



k 



l + A; 



-w^ + 



+ 



5!2/3„ir„^?"'^" ^°^ 
Ti + T2 + T3 + R^. 



3!2;g„ir„nl+27+Ai 7 

^T^^7 + ^7 
^^1+47 7 ' 



(4.13) 
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Proof of part (1): we consider 7 = 1/6, Ai 



with A = 



1 

7573 



and ip{x) 



KciM 



X — 



1/3 and A2 = 2/3 and get 



bC 



(4) 



■x^ — 



G 



(6) 



,(0) 



-X 



where we used Lemma 



exphcit formulas for ci 



4.6 



Ci(/3„,id7 C2 



Note that iE(H/,^^/(H/,/g)) 



and i?i/6 = -R/3„,K„,i/6 + 



cia:+C2X^+C3X^ with 
C4 



C2(/3„,i^„), C3 = csiPn^Kn) and C4 = 
we will compare the distribution 

of W^i/6 with a distribution with Lebesgue-probability density proportional to expf — ^^^^-^^ 



C4(/3„, K^, E{Wf/,)MWt/e)MW^/e))- Applying Theorem 



^^^^ — ^'^c ) • This density as well as the density p with ■?/' = p' /p fulfil assumption (B) , see P^ 
Lemma 2 .2] G iven the linear regression condition we are able to compute the bound give n in 
3.1 , Since \Wi/q — W[,q\ < ^^ we obtain with A 



Theorem 



n 



"^/^ and the help of Lemma 



3A 



\[ A ) ' 2 



E 



HWi/,) 



Exactly as in the proof of part (3) of Theorem 2.1 we have 

^2 



4.6 



1/6 



2A 



V 



E 



{Wy,~W'Y\Wy^ =0 n-1/3 



which completes the proof of part (1) where all regions A, B and C influence the limit distri- 
bution. 

Proof of part (2): We consider the linear regression condition 

k 1 



E[w^-wm 



K,{f3„ 



n 



l+A: 



-W^ + R^ 



with R^ = T2 + T3 + RJ, where R^ is defined in (4.13). Hence A 



1 



and ip{x) 



KcWn) 



X 



—x/cr and we compare the distribution of W^ with a N{Q,E{W^)) distribution. With \Wj — 



1 



K\<v^ 



and A :- 



1 



4.7 



^^x^ we have a'^l.5A^jE{W^) = 0{n^~'^) and A^/X = n^''^. From Lemma 
we know that |/?/3„,x„,-y| = 0{n^-'^) and therefore 0(^^) = r2-"^i'^(2-7,i+47,i+27+Ai)^ 



Lemma 4.6 With the proof of part (2) of Theorem 2.4 we have 

,2 



usmg 



a 
2A 



V 



E 



{W, 



w;y\w^ 



O (n-V2) 



and therefore the leading order is given by the order of A"^y/E[i?^], which leads to the three 
different cases, solving the minimization problem min(2— 7, 1+47, 1+27+A1) for 7 G (1/4, 1/2) 
and Ai > 0. 

Proof of part (3): We have exactly the same situation as in part (2). Therefore one has 
to solve the minimization problem min(2 — 7, 1 + 47, 1 + 27 + Ai) with 7 G (1/6, 1/4] and 
Ai > 2A2 — 1 = 1— 47. This leads to the four cases stated in the Theorem. 

Proof of part (4): With 7 = 1/6, Ai > 1/3 and A2 > 2/3 We consider the linear regression 
condition 



E[Wy^-W[,^\F_ 



1 

5!2/3„ir, 



G 



(6) 



(0) 



n 



5/3 



w!„ + ^1/6 



with _Ri/6 = Ti + T2 + -Ri/6, where Ri/^ is defined in (4.13). Hence A 

ip{x) = 



5/3 



and 



5!2A 



ri<"^^I,K„ (0) ^^ ^^d W6 compare the distribution of Wi/e with a distribution 
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with Lebesgue-density proportional to expf — gj^^^j. As in the proofs of part (2) and (3) we 

see that the leading order is the order of A^^i/E[i?^,g]. We now have to solve the minimization 

problem min(ll/6, Ai + 4/3, A2 + 1) for Ai, A2 > which leads to the result stated in the 
Theorem. 

Proof of part (5): With 47 = 1 — Ai, 7 G (1/6, 1/4) and 2A2 > Ai + 1 we consider the linear 
regression condition 



with R^ = Ti + T3 + i?^, where R^ is defined in (4.13). Hence A = 1+27+Ai and ip{x) = 
— 3,^i ^ ^isl,K„ (0) x'^ we compare the distribution of W^ with a distribution with Lebesgue- 
density proportional to expi — ^^^^^A . Again as in the proofs of part (2) and (3) we see that 



the leading order is the order of A ^i/E[i?^]. We now have to solve the minimization problem 
min(7, A2 - 1 + 27, 67 - 1) for 7 G (1/6, 1/4). 

Proofs of part (6) and (7): With 7 = 1/6 and Ai = 1/3, A2 > 2/3 or Ai > 1/3, A2 = 2/3 
we consider the the linear regression condition 

E|»Va - Wl,,m = ^G':'.. (0) J^W,% + ^/£.,<. (0) ^.n. + «V. 

with ^1/6 = Ti + Ri/Q or 

m.>. - wi„m = j^^^,w„, + ^g£U-.. (0) j^s^o + ^v« 

with i?i/6 = T2 + -Ri/6. In the first case solve min(ll/6, 1 + A2), in the second solve 
min(ll/6,Ai+4/3). 

Proof of part (8): Finally we consider the linear regression identity 
with R^ = T2, + R^. Solve min(7, 27, 67 — 1) = min(7, 67 — 1) for 7 G (1/6, 1/4). D 
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